Spin relaxation of conduction electrons in (llO)-grown quantum wells 
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The theory of spin relaxation of conduction electrons is developed for zinc-blende-type quantum 
wells grown on (llO)-oriented substrate. It is shown that, in asymmetric structures, the relaxation of 
electron spin initially oriented along the growth direction is characterized by two different lifetimes 
and leads to the appearance of an in-plane spin component. The magnitude and sign of the in-plane 
component are determined by the structure inversion asymmetry of the quantum well and can be 
tuned by the gate voltage. In an external magnetic field, the interplay of cyclotron motion of carriers 
and the Larmor precession of electron spin can result in a nonmonotonic dependence of the spin 
density on the magnetic field. 

PACS numbers: 72.25.Rb, 72.25.Fe, 73.21.Fg 



I. INTRODUCTION 



The long and tunable spin lifetime of carriers in semi- 
conductor low-dimensional structures is a crucial fac- 
tor for spintronic applications. Of particular interest 
in this context are quantum wells (QWs) grown from 
zinc-blende- type semiconductors on (llO)-oriented sub- 
strates. In such structures, the spin lifetime of conduc- 
tion electrons can be as long as several nanoseconds at 
room temperature"^ 3 - and tens of nanoseconds at low 
temperature 4 allowing for a long-range spin transport. 5 
Moreover, the spin lifetime can be widely tuned by the 
gate voltage 2 ^! 7 - and modifying the doping profiled These 
features are attributed to suppression of the D'yakonov- 
Perel' (DP) spin relaxation mechanism in symmetric 
(110) QWs. 9 The DP mechanism (Ref. pjj) usually lim- 
iting the electron spin lifetime is based on the precession 
of electron spins in an effective magnetic field induced by 
spin-orbit interaction in noncentrosymmetric media. In 
(llO)-grown QWs, the effective magnetic field caused by 
bulk inversion asymmetry (BIA) points along the growth 
direction.— Therefore, electron spins oriented along the 
QW normal do not precess in the field and the DP mecha- 
nism gets ineffective. If, however, the QW is asymmetric, 
structure inversion asymmetry (SIA) leads to the Rashba 
effective magnetic field^ which lies in the QW plane and 
speeds up the spin dephasing. Thus, by measuring the 
spin lifetime of carriers one can conclude on the Rashba 
field strength. 

Here, we develop a microscopic theory of the electron 
spin relaxation in low-symmetry two-dimensional struc- 
tures, such as QWs grown on (110)- and (113)-oriented 
substrates. We show that the mentioned above common 
analysis of spin dynamics can only be used to illustrate 
the increase of spin lifetime in symmetric (110) QWs. 
However, the simplified model misleads in describing the 
spin dephasing in asymmetric structures and determin- 
ing the Rashba field. In low-symmetry QWs, the growth 
direction is not a principle axis of the spin relaxation rate 
tensor. Therefore, the relaxation of electron spin initially 
oriented along the QW normal is described by a sum of 
exponential functions with different decay rates and leads 



to the appearance of an in- plane spin component. The 
magnitude and sign of the in-plane component arc deter- 
mined by structure inversion asymmetry of the QW and 
can be tuned by the gate voltage. We also study the ef- 
fect of an external magnetic field on spin relaxation and 
show that the interplay of the cyclotron motion of carri- 
ers and the Larmor precession results in a nonmonotonic 
dependence of the electron spin on the magnetic field. 

It is worth noting that the coupling of the in-plane 
and out-of-plane components of electron spin considered 
here is a feature of structures of low point-group symme- 
try and does not occur in (OOl)-grown (even asymmetric) 
QWs. In (OOl)-oriented structures, both the Rashba and 
Dresselhaus fields lie in the QW plane and their con- 
structive or destructive interference leads to the in-plane 
anisotropy of spin relaxation . 12 ' 13 ! 14 ' 15 By contrast, in 
(110) QWs the Rashba and Dresselhaus fields are orthog- 
onal and cannot compensate nor strengthen each other. 



II. SYMMETRY ANALYSIS 

Let us first consider the case of zero magnetic field. 
The time evolution of the spin density S(t), provided that 
the spin lifetime is longer than the carrier thermalization 
time, is described by 9 



dS a (t) 
dt 



(1) 



where G is the spin generate rate, e.g, due to optical 
pumping with circularly polarized light, T a p are com- 
ponents of the spin relaxation rate tensor T, and a, (3 
are the Cartesian coordinates. The form of the tensor T 
depends on the spin relaxation mechanism and the QW 
point-group symmetry. 

Asymmetric quantum wells grown on (HO)-oriented 
substrate as well as (113)-grown QWs are described by 
the point group C s which contains only two symmetry 
elements: identity and the mirror plane perpendicular 
to the QW plane 4^ In the case of (llO)-oriented struc- 
tures, the mirror plane is normal to the in-plane axis 
x || [110] and contains the axes y || [001] and z || [110]. 
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FIG. 1: (Color online) Arrangement of the principle axes x, 
y, and z of the spin relaxation rate tensor. Due to difference 
in the relaxation rates along the axes y and z, the relaxation 
of electron spin So \\ z leads to the appearance of the in-plane 
component S y (t). 



The symmetry analysis shows that non-zero components 
of the spin relaxation rate tensor in such QWs are T xx , 
^yyi ^zzi and T yz — ^zy- 

The presence of the off-diagonal components T yz and 
T zy indicates that both the in-plane axis y and the QW 
normal z are not principle axes of the spin relaxation rate 
tensor. Therefore, the decay of S y (as well as S z ) cannot 
be described by a single spin lifetime. The principle axes 
x, y, z and eigen values ji of the tensor T, i.e., relaxation 
rates along these axes, are found from the determinant 



det(T - 7 J) 



. 



(2) 



where J is the unit matrix 3x3. Solution of Eq. ^ has 
the form 



7* 
72 



(3) 



r^ + r 2Z + v /(r w -r 2Z )2 + 4rg z j /2, 

^-yy Tzz — \J (J^yy — Tzz) 2 + 4r^ z /2 , 



where the axes x, y, z are obtained from x, y, z by rotat- 
ing the coordinate frame around the axis x by an angle 
9, see Fig. Q] The angle 9 is given by 



tan 6 = 



2IY 



r 



yy 



(4) 



Solution of Eq. (fTJ) in the coordinate frame (x, y, z) 
can be then derived by projecting S{t) and G onto the 
axes x, y, and z. This calculation shows that, for pulse 
excitation, the time evolution of electron spin Sit) after 
the pulse has the form 



S x (t) = S Qx e-^\ 

S y (t) = S 0y (cos 2 9c-^ t +sin : 

+ S 0z cos 9 sin 9 (e" 7 ** - e" 7 ^) , 

S z (t) = S Qz (cos 2 9e" 1Bt +sin 2 6»e- 7it ; 

+ S 0y cos6»sm6»(e- 7i * - e" 7i *) , 



(5) 



where So is the spin density at t = 0. One can see that 
the relaxation of electron spin initially directed along the 
QW normal, i.e., when So || z, is described by two dif- 
ferent rates 72 and j y . Moreover, the relaxation leads 
to the rotation of spin in the (yz) plane resulting in a 
non-zero value of S y . Such a rotation is illustrated in 
Fig. [1] where the electron spin S(t) is shown as a sum of 
the components S y (t) and Sg{t) along the principle axes 
y and z. At t = 0, the spin is oriented along the QW 
normal, therefore, So y = So sin 9 and S02 — So cos 9. In 
the course of time, the components S y (t) and S z (t) de- 
cay at different rates, 75 and 75, respectively. Thus, at 
t > 0, the ratio Sy(t)/S z (t) is not equal to the initial 
one Soy/Soz- It means that the spin Sit) does not point 
along the QW normal anymore and has non-zero compo- 
nent S y (t) as shown in Fig. [TJ Similarly, the relaxation 
of electron spin oriented along the y axis leads to the 
appearance of out-of-plane component S z (t) . 

Besides experiments with time resolution, spin phe- 
nomena are widely studied in the regime of continuous- 
wave (cw) pumping, where the spin generation rate G 
is constant on the spin lifetime scale. The steady-state 
spin density in the regime of cw pumping can be directly 
obtained from Eq. (frj and has the form 

G ' (6) 



Sz = 



r 

1 XX 
G y T z 



GzT 



yz 



yy L zz 
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yy 
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yz 
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T yy T z 
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yz 



Equations ([5]) and ((6|) are general and describe the spin 
properties of (110) QWs no matter which microscopic 
mechanism determines the spin dephasing. 



III. D'YAKONOV-PEREL' MECHANISM 

In a wide range of temperature, carrier density and 
mobility, the spin lifetime in two-dimensional semicon- 
ductor structures is limited by the DP mechanism. 10 In 
this mechanism, components of the spin relaxation rate 
tensor in the collision-dominated regime are given by^ 

df(e k ) 



^afi — — 



Tl 



/(0) de k 



where t\ is the isotropization time of spin density, /(£&) 
is the distribution function of carriers, e k — Ti 2 k 2 / (2m*), 
to* is the effective mass, tt k is the Larmor frequency 
corresponding to the effective magnetic field caused by 
spin-orbit coupling, 5 a p is the Kronecker symbol, and 
the angle brackets denote averaging over the direction of 
the wave vector k. The Larmor frequency corresponding 
to fc-linear effective magnetic field in (110) QWs of the 
C s point group has the form 
2 

J~ifc = -riaiky, -a 2 k x ,(3kx) ■ (8) 
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FIG. 2: (Color online) Time evolution of the spin components 
Sz(t) and S y (t) calculated for different ratios ai//3 and the 
initial spin So oriented along the QW normal z. 



The parameter (3 originates from bulk inversion asymme- 
try while a\ and a 2 are non-zero only in QWs with struc- 
ture inversion asymmetry. We also note that a\ and a 2 
are linearly independent in QWs of the C s point group. 
The difference between a± and a 2 cannot be obtained in 
framework of the Rashba model which gives a\ — a 2 . To 
obtain the difference in a microscopic calculation of the 
band structure one needs to take into account both QW 
asymmetry and the lack of an inversion center in the host 
crystal^ 

It follows from Eqs. (J7J) and © that the components 
r ' a/ 3 in asymmetric (llO)-grown QWs assume the form 



T xx = (al + p 2 ) C , T yy = {a\ + p 2 ) C , 
T zz = {a\ + of) C , T yz = T zy = a 2 (3 C . 



(9) 



where the parameter C depends on the temperature and 
carrier density, and for the degenerate two-dimensional 
electron gas it is given by C = {Atiwi* /Ti A )Ep with Ep 
being the Fermi energy. Thus, for the DP mechanism we 
derive 7£ = {a 2 2 +(3 2 )C, j y = --{al +a 2 .+(3 2 )C : j z = a\C, 
and tan# = a 2 /(3 (see Ref. [18| for the case of ct\ = a 2 ) 
. One can see that all three spin relaxation rates are 
different and the rate 72 differs from T zz approximately 
by a factor of 2 even in the case of small angles 9. 

Figure [5] shows the time dependences of the spin com- 
ponents S z (t) and S y (t) given by Eq. for different 
ratios a±/p. It is assumed that the spin is initially ori- 
ented along the QW normal, Sq || z, and a± = a 2 . The 
time t is measured here in units of the in-plane relaxation 
time in symmetric QWs r sym = (3 2 C. The dependences 
plotted in logarithmic scale demonstrate that S z (t) is 
not described by one exponential function and S y can 
be comparable to S z even in the case of spin pumping 
along the QW normal. We also note that the sign of S y 
is determined by the sign of the product a 2 (3. Therefore, 



one can control the in-plane spin component by tuning 
the Rashba field, e.g., by the gate voltage. 



IV. SPIN RELAXATION IN MAGNETIC FIELD 

Now we consider the effect of an external magnetic field 
on spin relaxation. The effect is dual. Firstly, the mag- 
netic field causes the precession of nonequilibrium spin re- 
sulting in a spin depolarization (Hanle effect). Secondly, 
the magnetic field leads to the cyclotron motion of free 
carriers, which changes the direction of the wave vector 
k and thereby suppresses the DP spin relaxation mech- 
anisiri ] 19 ' 20 i 21 i 22 Both effects can be fruitfully treated in 
the framework of the spin-density-matrix formalism. In 
this approach, the electron distribution is described by 
the spin density matrix Pk = fk + Sk ■ cr, where fk and 
Sfc are the functions of particle and spin distributions in 
the k space, respectively, and a is the vector of the Pauli 
matrices. The spin distribution Sk, within the relaxation 
time approximation, satisfies the kinetic equation (see, 
e.g., Refs. [HIH) 



-^-+Skx(n L - 

ot 



ds k 

dk a 



k x 



eB 



Sk - (Sfc) 



( 10 ) 

where $7^ is the Larmor frequency whose components are 
given by flL.a = PBQapBp /Ti, hb is the Borh magneton, 
g a /3 is the effective g-factor tensor in the absence of cy- 
clotron motion, e is the electron charge, c is the speed of 
light, and Gk is the spin generation rate into the state 
k. We note the g- factor can be strongly anisotropic in 
quantum well structuresj 23 ' 24 i 25 In particular, in (110)- 
grown QWs of the C s symmetry the components g xx , 
gyy, gzz, gyz, and g zy are linearly independent and can 
be non-zero. 

The anisotropic part of the spin distribution function 
5sk = Sfc — (sfe) is caused by spin-obit splitting of electron 
states. In the collision-dominated regime, i.e., at fifcTi <C 
1, it is much less than the isotropic part (sk) and, to the 
first order in OfcTi, has the form 



Ss k 



ti [(sk) x (Ofc + m c Tid£lk/d(pk)] 

1 + (^ c Tl) 2 



(11) 



Here, lo c = eB z /(m*c) is the cyclotron frequency, ipk is 
the polar angle of the wave vector fc, ipk = arctan(fc y /fc a: ), 
and (sfc) is the quasiequilibrium spin distribution which, 
in the case of small degree of spin polarization of a two- 
dimensional electron gas, is given by2> 



(s k ) = - 



2nh 2 df(ek) 
m*/(0) de k 



S . 



(12) 



Equation (fTTj) is derived taking into account the fact that 
Jlfc is linear in k and assuming that Qlt~i <C 1 and Gk 
is an isotropic function of k. 

Finally, substituting expressions given by (fTTj) and flj 
for Ssk and (s/c), respectively, and summing up Eq. |TI 
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FIG. 3: (Color online) Dependences of the electron spin 
S z on magnetic field aligned along the QW normal for dif- 
ferent ratios ot\/j3. Solid and dashed curves correspond to 
^L,zT sym /LJ c Ti = 5 and r2L jZ r sym /aj c ri = 1, respectively. 



over the wave vector k, one obtains the kinetic equation 
for the total spin density 



S a (t) 
dt 



[S(t)xn' L ] a =G a -J2KpSp(t), (13) 



where T' a g — T a p/[1 + (w c ri) 2 ] are components of the 



spin relaxation rate tensor in the magnetic field, fl' L = 
Ct L + SCIl, and 5tt L is a contribution to the precession 
frequency caused by cyclotron motion of free carriers in 
QWs with spin-orbit splitting (see Ref. [19]), 



sn, 



LO c Tl 



1 + (w c Tl) 2 



n df(e k ) 
2/(0) de k 



n k x 



dfl k 



de k 
(14) 

In fact, this contribution can be considered as a correc- 
tion to the effective g-factor which depends on the relax- 
ation time T\ and the magnetic field strength. We note 
that the frequency SSIl is expressed via the same pa- 
rameters which determine the spin relaxation rate tensor 
and is comparable to T' a a in the field lo c t\ ~ 1. Taking 
into account the form of £l k given by Eq. ((8]) , we obtain 
non-zero components of 5£Il in (llO)-grown QWs 



l + {LU c Tiy 



1 + (w c Ti) 2 



C. (15) 



Equation (|13[) describes the spin dynamics in external 
magnetic field. It yields that the dependence of the out- 
of-plane spin component S z , which is usually studied in 
experiment, on the magnetic field in the case of cw spin 
pumping along the QW normal has the form 



S z = 



(v r 

V xx w 



x XX 



(V V 

y yy * 



yz) 



12 



1 yy lL L, v 



_i_ on o' o' 



(16) 



Shown in Fig. [3] are the dependences of S z on the field 
B || z (Faraday geometry) plotted for different ratios 
ai//3 and 0, LtZ T sym /u c Tx, «i = "2 and g yz = 0. One 
can see that the magnetic field dependence of S z can 
be nonmonotonic even in the Faraday geometry and the 
effect is more pronounced at |nL iZ r sym /a; c ri | ^> 1. Such 
an unusual behavior is attributed to the noncoincidence 
of the growth direction with a principle axis of the spin 
relaxation rate tensor, see Fig. [TJ As it was discussed 
above, in zero magnetic field, the electron spin initially 
oriented along z is decomposed into the projections S z 
and Sy which decay at different rates. Since 7^ > 7^ 
for the DP mechanism, the component Sy rapidly decays 
and the spin polarization is determined by the long-life 
component S^. The external magnetic field deflects the 
electron spin from the slow-relaxation axis z speeding up 
the spin dephasing. Thus, the Larmor precession leads to 
the decrease of spin polarization in the magnetic field. At 
stronger fields, the cyclotron motion suppresses the DP 
mechanism of spin relaxation, which results in a growth 



of S z . The interplay of the Larmor precession and the 
cyclotron motion of free electrons in asymmetric (110) 
QWs leads to the nonmonotonic dependence of the spin 
density S z on magnetic field as shown in Fig. [3] 

In summary, it is demonstrated that the in-plane and 
out-of-plane components of electron spin are coupled in 
quantum wells of low space symmetry. This opens up 
additional opportunity to manipulate electron spins in 
low-dimensional semiconductor structures. 
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